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Abstract 

The BCR-ABL translocation is found in chronic myeloid leukemia (CIVIL) and in Ph+ acute lymphoblastic leukemia (ALL) 
patients. Although imatinib and its analogues have been used as front-line therapy to target this mutation and control the 
disease for over a decade, resistance to the therapy is still observed and most patients are not cured but need to continue 
the therapy indefinitely. It is therefore of great importance to find new therapies, possibly as drug combinations, which can 
overcome drug resistance. In this study, we identified eleven candidate anti-leukemic drugs that might be combined with 
imatinib, using three approaches: a kinase inhibitor library screen, a gene expression correlation analysis, and literature 
analysis. We then used an experimental search algorithm to efficiently explore the large space of possible drug and dose 
combinations and identified drug combinations that selectively kill a BCR-ABL+ leukemic cell line (K562) over a normal 
fibroblast cell line (HVlR-90). Only six iterations of the algorithm were needed to identify very selective drug combinations. 
The efficacy of the top forty-nine combinations was further confirmed using Ph-I- and Ph- ALL patient cells, including 
imatinib-resistant cells. Collectively, the drug combinations and methods we describe might be a first step towards more 
effective interventions for leukemia patients, especially those with the BCR-ABL translocation. 
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Introduction 

The BCR-ABL fusion protein results in the Philadelphia 
chromosome and is present in 95 % of chronic myeloid leukemia 
(CML) patients and 20-40% of adult acute lymphoblastic 
leukemia (ALL) patients. [1-3] The BCR-ABL translocation is 
characterized by a constitutively active fusion tyrosine kinase and 
plays a causal role in CML. It is therefore an attractive target for 
drug therapy. [4,5] Imatinib mesylate (Gleevec/ Glivec), a selective 
Abl kinase inhibitor, has been shown to have marked efficacy as a 
single agent and is used as a first-line therapy for the treatment of 
CML and Ph-H ALL patients. [6] 

Despite the remarkable hematologic and cytogenetic responses 
to imatinib, a substantial number of patients develop resistance, 
especially in advanced cases, leading to disease relapse. [7] 
Resistance to imatinib has been attributed to secondary genetic 
mutations in the target and to downstream mechanisms. [8] In 
addition, imatinib is unable to eradicate the leukemic blast cells, 
thus only inducing remission. [7] To circumvent these problems, 
many research groups have investigated combinatorial drug 
regimens with drugs that have different modes of actions. [9-1 1] 
Most of the combinatorial studies, however, considered combina- 
tions of two or three drugs, which may not suffice to achieve 
complete inhibitory action on multiple oncogenic pathways. [12] 



Results 

Search flow for optimization of combinatorial drugs 

A scheme of our search strategy to identify optimal drug 
combinations is illustrated in Figure 1 . We first selected 1 1 drugs to 
be combined with imatinib and then used the BCR-ABL-I- K562 
cell line as a model system. The drug combinations were then 
optimized using iterative loops that process input data from 
experimental readouts. We applied an algorithm suggesting 
combinations to be tested at each iteration. 

Dose responses of single agents were performed to determine 
the starting doses to be used. We then performed a pair analysis, in 
which all the possible pair combinations of agents were studied; 
indicated in the Figure 1 as "triplets" because imatinib was also 
added, at fixed dose. The data from the pairH-imatinib analysis 
(triplets) were fed into a model that estimates the effectiveness of 
large drug combinations, assuming only pairwise interactions 
among compounds. The top combinations predicted by this model 
were then used in the iterative loop as initial conditions for the 
experimental search. At each iteration new combinations were 
generated and experimentally tested, and this was repeated until 
effective drug combinations were identified. The best combina- 
tions identified using the K562 cell line were vahdated using 
primary cells from BCR-ABL-I- ALL patient xenografts and ALL 
patients. 
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Figure 1. Scheme of our search strategy (a) anti-leukemic drugs were selected using three approaches: kinase inhibitor library screens, 
correlation analysis, and literature survey, (b) Dose responses of single agent and pair-wise analysis with a fixed dose of imatinib (pairs and triplets) 
were performed, (c) The drug combinations were optimized using iterative search algorithm (d) The optimal drug combinations were validated with 
primary patient cells. 
doi:10.1371/journal.pone.01 02221. gOOl 



Selection of anti-leukemic drugs 

We employed three difTerent methods to select the set of 1 1 anti- 
leukemic drugs that are studied with imatinib in the combinatorial 
drug searches. We first screened an EMD kinase inhibitor library 
that consists of 244 kinase inhibitors with/without imatinib 
(Figure 2). The cell viabilities were measured by ATPlite assays 
in the leukemic cell line (K562) and in a control fibroblast cell line 
(IMR-90) [13,14] 72 hours after the treatment with the kinase 
inhibitors. We have rank-ordered the kinase inhibitors according 
to the selectivity index, which is defined as the ratio v„/V;., where 
Vn and Vj- indicate the viability of IMR-90 and K562. Four highly 
ranked kinase inhibitors were selected based on the analysis; A15, 
CI 6, E3, 122 (Table SI). We also selected three kinase inhibitors 
that showed synergistic or additive effects when tested with 
0.125 |J,M imatinib based on the imatinib combination index (v^./ 
Vmiatinib) whcrc v^ and Vin,atmii) indicate the viability of K562 
without/with imatinib); A07, A12, P15 (Table SI). This fixed dose 
of imatinib was chosen so that it gives only 10-20% cell killing 
according to the dose response curve of imatinib. 

In addition, we chose ABT263, Axitinib, and 17AAG based on 
a correlation analysis between drugs and imatinib responses. The 
correlation signs were derived from the Pearson correlation 



between mRNA expression and drug IC-,o values on cancer cells 
data from Garnett et al. [15] The imatinib- 1 7AAG pair was among 
the top scoring negative-positive and positive-negative pair set, 
whereas ABT263 and Axitinib were among the top of the positive- 
positive and negative-negative pair sets (Table S2). The rationale 
for this analysis was to identify drugs that might either share gene 
expression correlations (the positive-positive and negative-negative 
group) with imatinib or that might have opposite correlations (the 
negative-positive and positive-negative set). We hypothesized that 
the first type of compounds might act on similar pathways, via 
different targets, and therefore directly reinforce the effects of 
imatinib, while the second type might ideally complement it. 
These strategies should counteract resistance mechanisms of 
diflFerent types. 

We also added the mitogen-activated protein kinase (MEK) 
inhibitor PD0325901 to the drug list based on literature data. 
[16,17] The list of the small molecules used for our combinatorial 
study and their targets are shown in Table 1. 

Dose responses and pair analysis 

Dose responses for the selected small molecules were measured 
over a 100-fold concentration range (seven-points, as shown in 



PLOS ONE I www.plosone.org 



2 



July 2014 I Volume 9 | Issue 7 | e102221 



Drug Combinations Containing Imatinib for BCR-ABL+ Leukemias 



(a) 



(b) 



(c) 



' 1163840 
1150400 
1139CBQ 
1051580 
1OG5660 
427060 
1024960 
1094120 



999S50 U63140 1045340 342480 

1059200 1084360 1106740 10M520 

1036320 1115480 864340 993450 

993140 59960 1043120 9904S0 

965540 | Sliiol 934340 967680 



105322o | 105420| 

10S8380 1075920 

1035140 942520 

1026000 905920 

897400 914320 

909920 920540 



962160 
979600 
941780 
945620 



902280 953960 945040 

974400 978560 967600 

8^860 94^20 1046300 933640 

820620 961080 7864S0 1004820 

989540 1054140 93800D 

102O1BO 995920 975140 

1049940 965020 947340 

1135880 1D34400 1053S4O 



920320 
918520 
917180 
9334S0 



lOO&WO 875380 
861460 S583C0 



1047360 
1035320 
1057020 
[ 1133130 
^ 1098900 



63220 824980 985480 858880 

792560 1003980 822400 

980740 921140 861620 

533S60 999100 845920 



lllSOQO 1037440 
1100010 958920 



991920 1037860 



973120 917120 
1099620 920520 



875340 
336160 
901280 
S7310O 
166730 
924040 
375240 
933980 
947460 
992780 



9 10 

349560 1CO49O0 

931000 994560 

350520 1C02380 

9083SO 953730 

452960 874420 

9D9280 967760 

845460 846340 

347360 910760 

755380 60700 

865200 904560 

329250 950800 

351880 913960 

911750 800200 

925100 1040540 

84360 963040 

924920 1031580 



S0 | 314460| 



14 



IB 



952640 1009280 

953000 354430 

914240 910400 

360740 843600 

550080 922160 

887460 830930 

853200 873440 

50240 899520 



81100 

853620 
396400 



9101Bq | 293eo| 

959040 71580 

983860 74752cj^ 

931030 9034D0 

883300 343060 

870540 790320 

822200 751620 

897440 823200 

767530 664320 



21 



22 



1004300 975450 126540 390240 
988160 1005000 852080 953380 
925680 992580 90^60 



853130 
851760 
804740 
818920 



818940 
753480 
867000 



977220 
950860 



952160 
931330 945540 



935620 1016560 870200 954500l 



811440 
771040 



) | 80768C| 



101140 

875520 
16D7DO 
958240 



336S0O 
337500 



873360 
912780 
831940 
820150 
933320 



912420 
770760 
S73320 
763580 
858180 
347380 
908220 
939320 
905500 
837850 
951540 
W8160 



900920 
927460 



437360 
872560 
835330 
326340 
301340 
901800 
848680 
133440 
822900 
945120 
520550 
577520 



966280 9225C0 

911060 987680 

953500 975420 

965040 1023680 

954260 951460 

906700 920320 

859520 910080 

920180 904960 



356640 
864740 
904740 

917320 



879580 866460 



23 
1023540 
1004150 
1D12540 
1011560 
953240 
959620 
984000 
946500 
915630 
919430 



W65C0 958030 

890360 855780 883980 1005140 

930880 ^6800 972350 990420 

130900 908160 1002900 

924160 976900 9S3130 



695220 
728&40 





1 


2 3 


4 5 


6 7 
85154o| 48560| 
948200 978920 


8 


9 


10 


11 12 13 
753140| 1357B0| 345280 
843200 877620 72160 


A 
B 


935800 
904740 


105Z6O0 952100 
1102460 869020 


980160 328500 
976400 852720 


903120 
924100 


676820 
804000 


871200 
900160 


C 
D 

j E 


930140 

914040 
347950 


1081150 945740 
106S620 381020 

io6592o|~ Agesoj 

1094160 91^20 


847140 S78520 
910260 928CO0 
394030 3^960 
918420 S87S40 


869660 855380 
912960 834640 
905960 835350 
897600 821240 


869120 
867000 
526320 
310760 


793820 
848000 
454560 
814460 


874780 
829640 
S47520 


833400 340080 791680 
827200 891460 830600 
762150 789120 793420 
519260 845750 735700 


G 


870540 
824140 
743480 
549320 


1096260 851900 
1007720 855040 
1073300 903740 
1066140 817200 


907680 913700 
910400 895460 
881200 870400 
917520 929800 


855040 737300 
890140 745000 
758460 828920 
837420 810720 


306940 
769620 
301420 
8044S0 


735680 
305540 
666000 
757250 


746460 
779330 
4E180 

777100 


803000 815680 662240 
790820 859400 667840 
42020 779080 529960 
754800 840000 647380 


K 
M 


816630 
873120 
SS4960 
894930 


1105200 9004CO 
1143720 912400 
112^00 923140 
1145360 978140 


57640 803440 
901730 570020 
847520 900430 
975360 409820 


913920 792080 
942300 774340 
939660 80560a 
931060 793780 


215740 

327040 
361140 
919120 


796740 
752160 
849180 
352860 


766300 
822340 
773940 
915420 


51C840 714630 760040 
853S20 845950 594540 
795420 779380 751300 
730260 905330 303250 


0 
P 


934620 
920140 


1176420 979220 
1066420 947020 


890820 1020880 
1046560 1017440 


1C01020 917320 
1C»5300 831220 


952120 
970940 


nA2ao 

862020 


858380 
843200 


92SW0 874240 83SS40 
834420 913280 795340 



788040 676560 

628280 709180 

784240 725360 

648900 725520 



852100 
S30I60 



8554a) 798930 339440 
327030 825730 776610 
671220 786440 766220 
642160 767200 722060 
749440 
755620 
755010 
206340 
310140 



£0 374240 
00 870540 
S0 | 43738q| 



753M0 
824840 
115920 

8700)0 



18 

900920 
872640 
841440 
810060 
712050 
764920 
607700 
794940 
7776t» 



S230S0 
813520 
925040 



851520 
831150 1003040 
339200 945420 



19 
105100 
375440 
33^00 
334650 
409760 
739950 
305140 
304200 
763780 
314220 
761920 



372330 943520 
457550 1036340 
884160 961380 



S%960 943500 
940340 SS74SO 
SSS380 892100 
864720 914400 
493S0 873920 
862960 907900 
809940 772SOO 
835450 858720_ 
724020 S34t 
310180 841420 
S53CO0 834340 
862SS0 708420 
884700 873100 
827H)0 
193000 
927400 



22 
757920 
870OCO 
374140 
363360 
789240 
7500C0 
763140 
792320 



913140 
918480 



9B600 
947280 

931200 
933120 



905340 
900960 
334520 
930720 
939140 
935580 
908220 



925620 
930060 



907520 
861260 
814420 



7592C0 
9055C0 
802720 



810140 
853240 
845750 
890720 
950940 
946620 
918100 
895240 



337240 
336720 
862900 
362200 
853180 
828220 



IMR 

K562 

K562 with Imatinib 




1100000 
1000000 
900000 
800000 
700000 
600000 
500000 
400000 
300000 
200000 
100000 
SOOOOi 
250Qo! 

looooj 



RLU 



Figure 2. Selection of anti-ieui<emic drugs from icinase inhibitor library screens (a) A representative screening plate of a kinase inhibitor 
library consisting of 244 inhibitors. The ATP content of K562 cells (a measure of cell viability) was measured and the luminescence intensity of each 
well was color-coded according to the scale presented (RLU, relative luminescence unit). The corresponding compound ID in the plate map is shown 
in Table SI. (b) A representative screening plate of a kinase inhibitor library treated with 0.125 |iM imatinib. The kinase inhibitors selected from the 
screens are indicated with black borders, (c) The kinase inhibitors are rank-ordered from the lowest to the highest cell viability in K562 cells without 
imatinib (red bars). The corresponding viability of IIVIR and K562 cells with imatinib are indicated using blue and green bars. Several drugs are 
selectively killing K562 cells. 
doi:1 0.1 371 /journal.pone.01 02221 .g002 



Figure SI). The IC5Q values for most of the drugs ranged between 
0.1 and I |J,M although 122 and 17AAG showed 1 0-fold lower 
values of IC50 (Table S3). Additional dose responses in combina- 
tion with imatinib confirmed the additive effects of the three drugs 
we chose from the initial screening (A7, A12, PI 5). We also 
detected an additive effect of A15 at low doses. ABT263, selected 
from the correlation analysis, also exhibited additive effects. For 
the subsequent pair analysis, the doses that cause a K562 cell 
killing between 20—30% in the presence of imatinib were chosen as 
high doses whereas those with killing between 10-20% were 
defined as low doses. These doses were selected in order to avoid 
that the cytotoxicity of single agents would dominate the response 
to the combinations. The full triplet analysis, which comprises 242 
drug pairs with a fixed dose of imatinib, was performed and the 
highest selectivity achieved in these combinations was 2.8 for the 
combination of A15 and 17AAG plus imatinib (Table S4). 



Identification of optimal drug combinations using a 
search algorithm 

We have developed a search algorithm that integrates two 
stochastic combinatorial optimization methods known as particle 
swarm optimization and genetic algorithm [18]. Initially, five 
doses of each single agent were considered and coded as number 
'0' to '4' (Table 2). The level 2 and 3 correspond to the low and 
high doses in the pair analysis, respectively. During algorithmic 
searches we realized that the doses used for some of the drugs were 
too high, resulting in combinations that were too toxic to the IMR- 
90 control cells. For other drugs the dose of maximal selectivity 
was not clear. Therefore after the third iteration, the concentration 
grid was reshaped to include two higher or lower levels depending 
on the drug performances (Table 2). After this reshaping of the 
drug concentration grid the algorithm was able to find highly 
selective combinations. 

The initial iteration was based on the 242 pair-wise measure- 
ments with imatinib (triplets) and a geometric average model that 
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Table 1. The list and targets of the drugs used In the combinatorial studies. 



Drug name 


Well ID 


Pubchem ID 


Drug Targets 


Target family 


Effector pathway/biological process 


Akt Inhibitor IV 


A07 


5719375 


AKT 


AGC 


Apoptosis, Metabolism, PI3K/MT0R 


Alsterpaullone, 2- 
Cyanoethyl 


A12 


16760286 


CDKl, GSK3B, CDK5 


CMGC 


Mitosis 


PDKl/Akt/FIt Dual 
Pathway Inhibitor 


A15 


5113385 


PDK1/AKT/FLT 


AGC, ATYPICAL, TK 


PDKl/AKT/FLT 


Cdkl/2 Inhibitor III 


C16 


5330812 


CDKl/2 


CIMGC 


Cell Cycle 


EGFR Inhibitor 


E03 


9549299 


EGFR 


^^^^^ 


ERK Signalling, PI3K/MT0R 


JNK Inhibitor IX 


122 


16760525 


JNK 


CIMGC 


Stress Pathways 


WHI-P180, 
Hydrochloride 


P15 


5687 


CDK2 


CMGC 


Cell Cycle 


ABT263 


N/A 


24978538 


BCL2, BCL-XL, BCL-W 


Other 


Apoptosis 


Axitinib 


N/A 


6450551 


PDGFR, KIT, VEGFR 


RTK 


Angiogenesis, ERK Signalling, PI3K/MT0R 


17AAG 


N/A 


6505803 


HSP90 


Other 


Other 


PD0325901 


N/A 


9826528 


MEK 


STE 


ERK Signalling 


Imatinib 


N/A 


5291 


ABL, KIT, PDGFR 


CTK, RTK 


Cytoskeleton, ERK Signalling, PI3K/MT0R 
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can estimate the selectivity of an arbitrary combination (see 
Methods). The searcli algorithm used these estimates to suggest 
the initial 270 combinations that were experimentally tested. After 
this initial iteration, the best selectivity was already improved by a 
factor of four compared to the best triplet selectivity of the 
previous step. (Table S4). At the sbcth iteration, we had 
cumulatively tested 1607 combinations and the iterative search 
converged to an optimal combination set. 

Figure 3 illustrates the optimization of drug combinations. The 
function being maximized is selectivity. Note that after the third 
iteration, when the dose levels were adjusted, the selectivity 
drastically improved. In the last three iterations we also refined the 
selectivity measure using a "corrected selectivity" that further 
penalized combinations that are toxic to the IMR-90 control. The 
corrected selectivity was obtained through a multiplicative term 

(^l+exp[— , where Vo is a cut off value describing the 

maximum allowed toxicity on normal cells, vn is the viability of 
the normal cell, and a is a smoothing parameter. We were able to 
attain a selectivity of 145 in the sixth iteration and a few 
combinations showed viabilities of more than 90% for IMR-90 
control cells and of less than 2% for K562 cells (Table S4). The 
best combinations had relatively high doses of A7 and Axitinib and 
low doses of 17AAG and ABT263. 

After the final iteration, we selected the top 200 combinations 
and performed a cluster analysis based on the euclidean distance 
metric. We observed two clusters; a larger cluster of 146 
combinations and a smaller cluster of 54 combinations 
(Figure 3(c)). The major difference between two clusters is in the 
dose of Axitinib. The analysis of the top 200 optimized 
combinations also suggested the optimal dose ranges of the 
individual drugs. Most of the drugs were used at a concentration 
lower than 100 nM, which represents a significant reduction 
compared to the concentration range at which the single drugs 
were effective. For example, the highest selectivity we could 
achieve using single drugs was 12 with 0.05 |J.M 17AAG whereas 
the best combinations required a 1 0-fold lower dose of the same 
drug to reach 10-fold higher selectivity (Table S4). 



Confirmation of optimal drug combinations using an 
imatininb-resistant cell line and imatininb-resistant and 
sensitive primary patient cells 

To demonstrate that our optimal combinations can also be 
apphed to imatinib-resistant cells, we selected 49 drug combina- 
tions and tested them in the SUP-B15 cell line and ALL patient 
cells with and without the BCR-ABL translocation (Table S5). The 
combinations were selected using two criteria: assuring that the 49 
combinations were representative of both cluster 1 and cluster 2 
using a 70% viability of IMR-90 as cut-off and assuring that low 
toxicity combinations with 90% viability of IMR-90 as cut-oflF 
were included. Effective cell killing was also observed with the 
imatininb-resistant cell line, SUP-B15, which showed even higher 
selectivity than K562, assuring that drug combinations work well 
on both imatinib-sensitive and resistant cells. We were able to 
achieve a selectivity of 19 (BOOl) and 24 (B002) with selected drug 
combinations when tested in BCR-ABL+ ALL xenograft derived 
patient cells (Figure 4 and Table S6). In addition, two BCR-ABL+ 
patient-derived cells with the T315I mutation (BOOS and B004) 
were included in order to determine if the combinations are also 
effective in known imatinib-resistant patient cells. Interestingly, 
higher selectivity was achieved with one of the imatinib-resistant 
patient cells (B003) but the other patient cells also responded well. 
The data shown in Figure 4 also demonstrate that several drugs 
and combinations we identified with K562 cell line are effective 
against BCR-ABL- ALL cells although the selectivity was not as 
high as for BCR-ABL+ ALL cells or cell lines. 

A statistical analysis of the drug combination results shows that 
the mean differences between BCR-ABL+ ALL patients and 
normal subject are larger than those between ALL or Burkitt's 
lymphoma patients and normal subject (Table 3). The differences 
were statistically significant for all the BCR-ABL+ patient cells and 
for all but one of the BCR-ABL- patient cells. On the other hand, 
the single agents showed mosdy negative mean differences, 
indicating that the single agents/imatinib pairs were generally 
not effective at the doses used for the combinations. Additionally, 
CD 34+ hematopoietic stem cells (HSCs) were used to further 
study the toxicity of the drugs. We also observed high selectivities 
of BCR-ABL+ ALL ceUs against CD34+ HSC controls, although 
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the viability of CD34+ cells tend to be lower than IMR-90 cells 
(Table S7). 

ABT263 was, however, highly selective even as a single agent rn 
leukemic cells although combinations always performed better 
than individual drugs (Table S6). Interestingly, the patient cells 
responded differendy to some of the drugs. For example, only one 
patient (B002) responded to PD0325901 whereas the other 
patients did not respond to it, at least as a single agent. Primary 
cells isolated from Burkitt's lymphoma did not respond to the drug 
combinations as well as leukemic cells. 

Discussion 

In this study, we have identified a list of potential drugs for 
BCR-ABL+ leukemias using three different approaches (kinase 
inhibitor screens, correlation analysis, and literature survey) and 
optimized the drug combinations using an iterative experimental 
search algorithm. Among the identified drugs, some showed 
additive effects when paired with the BCR-ABL+ targeting drug 
imatinib. Within only six iterations, we were able to find effective 
drug combinations that induced very selective killing of the 
leukemic cells (K562) versus the control fibroblast cells (IMR-90). 
The selectivities we achieved with the more effective drug 
combinations were significantly higher than those with the single 
agents, both in cell lines and in patient cells. 

17-AAG is already known for its cytotoxic effects on imatinib- 
resistant patient cells and for its synergistic activities with imatinib 
in AML patient cells, in agreement with our correlation based 
analysis. [10,19,20] However, the doses of 17-AAG used in the 
previous combination study were much higher than those in the 
present study. [10] Some of the kinases we targeted in our 
combinatorial experiments, such as AKT and JNK, are known to 
be downstream of the BCR-ABL fusion protein. [21,22] The 
BCL2 family inhibitor ABT-263 was highly potent in our study, 
selectively killing both K562 cells and leukemic patient cells at nM 
ranges. A previous study reported that ABT-737, an analogue of 
ABT-263, greatly enhances the apoptotic effect of INNO-406, a 
second-generation BCR-ABL inhibitor, against BCR-ABL positive 
leukemic cells. [23] In that study 17-AAG showed synergistic 
effects with ABT-263. The use of these two small molecules was 
suggested by our computational correlation analysis, but they or 
their analogues are already under clinical evaluation for the 
treatment of leukemia. [24,25] In addition, a previously published 
report shows that subutoclax, a pan-BCL2 inhibitor, sensitizes 
leukemia stem cells to dasatinib treatment. [26] These literature 
results validate our methodology for selecting the set of anti- 
leukemic drugs to be used in the experimental searches. 

The best drug combinations selected through the search with 
the K562 cell line also showed great efficacy in patient cells. The 
drug combinations were especially effective on BCR-ABL+ 
leukemic cells, but they also appeared to have less pronounced 
but stiU strong cytotoxic effects on several BCR-ABL- ALL patient 
cells (Figure 4). This is not unexpected, because they target 
proteins generally implicated in leukemia proliferation, part of 
common signaling cascades between BCR-ABL+ and BCR-ABL- 
leukemic cells. For example, it was previously reported that ABT- 
263 was not only able to induce cell death in BCR-ABL+ CML 
but also to induce complete remission of BCR-ABL- ALL cells in a 
mouse model. [24,27] Although imatinib was used to target BCR- 
ABL, it has been reported that this compound also inhibits ARG, 
PDGFRA, PDGFRB, and c-KIT. [28,29] AdditionaUy, kinase 
profiling studies performed by Anastassiadis et al. show that many 
other kinases are targeted by imatinib. [30] Further studies will be 
required to understand in detail the mechanism of action of 
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Figure 3. Iterative searcKi for a highly selective drug combinations (a) Increased selectivity after each Iteration step. The lines represent the 
selectivity of the drug combinations at as triplets and at each of the six iteration steps, (b) A dot plot showing corrected selectivity (see Results) 
through the iterative searches. The pie symbols shows the relative contributions of the 1 2 drugs for the best combination obtained at each iteration, 
(c) Cluster analysis on the optimized drug combinations. The top ranked 200 combinations clustered into two groups and the frequencies of the 
doses of the individual drugs are shown for each cluster. 
doi:1 0.1 371 /journal.pone.01 02221 .g003 



Imatinib in BCR-ABL- cells. In addition, the drag combinations 
also proved highly effective in the imatinib-resistant cell line, SUP- 
B15, [31] and in BCR-ABL+ leukemic cells with T315I mutation, 
pointing to the potential application of these drug combinations in 
tyrosine kinase inhibitor-resistant patients. 

Our study presents a novel search algorithm for the identifica- 
tion of drug combinations that markedly increase the selective cell 
killing of leukemic cells. The highest selectivity achieved through 
iterative searches was not obtainable with a single agent or 
combination of two or three agents but required combining a 
larger number of drugs (see Table S4 and Figure 3). Whereas the 
current methodology to identify effective drug combinations often 
involves exhaustive testing, several papers have previously 
suggested approaches to combinatorial control using drugs based 
on biological search algorithms. [32,33]. Our group has previously 
shown that algorithms derived from the field of digital commu- 
nication can be used to find combinatorial therapies. [33] . Wong et 
al. developed a closed loop control algorithm using a microfluidic 
platform to implement an iterative stochastic search for optimal 
drug combinations. [32,34,35] In this study, we integrated particle 
swarm optimization and a genetic algorithm [18] to perform a 
biological search for drug combinations. We provide a detailed 



description and pseudocode in the methods section. Both methods 
are inspired by biological phenomena and have been very 
successful in multiple computational and engineering applications. 
[18] 

The particle swarm algorithm mimics the behavior of a colony 
of insects that explore a field for food: each insect locally explores a 
section of the field, but there is exchange of information between 
them. Therefore, if an insect finds a region with abundant food, it 
will "inform" other insects and will attract them to that region. In 
our case, we consider a set of "agents" moving in the multi- 
dimensional space of all possible drugs and doses. Each agent 
explores a "neighborhood" of combinations, defined by adding 
and subtracting one dose. When the algorithm evolves, each agent 
has two options that are randomly chosen: (1) it can remain within 
its neighborhood and move towards the direction of the local 
optimal combination, or (2) it can look at how successful have been 
other agents and move towards neighborhoods of "more 
successful" agents. 

The genetic algorithm is inspired by the concepts of crossover 
and mutation in natural selection. We can think of a "population" 
of drug combinations in which the genotype of each individual is 
represented by the sequence of doses corresponding to an ordered 
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Figure 4. Test of optimum drug combinations in the patient cells. The best 49 drug combinations or single agents used for the combination 
were tested in the cells from BCR-ABL+ and BCR-ABL- ALL patients. The selectivity (vs control IMR fibroblasts) obtained in each patient specimen was 
color-coded from red (lowest) to green (highest). The drug combinations were also tested in Burkitt's lymphoma patient and in normal subject white 
blood cells as controls. 
doi:1 0.1 371/journal.pone.01 02221 .g004 



set of drugs. Top individuals are allowed to "mate" by swapping 
50% of their dose sequence to generate new "children" 
combinations. A small probability of mutation is also added, 
meaning that the dose of a randomly chosen drug is modified. In 
our procedure, we first identified combinations according to the 
particle swarm method, and then, using the top 10% individuals, 
we generated "children" combinations until a maximum value 
determined by the number of plates available is reached. 

The most effective combinations were composed of most of the 
drugs at low doses. These large combinations are analogous to the 
combinatorial strategies used by nature, as shown by the large 
number of regulatory molecules (e.g. transcription factors, kinases 
and micro RNAs) that determine cell fate and differentiation. [12] 
Large drug combinations are also more likely to overcome or 
prevent cancer resistance, which is one of the main factor limiting 
the efficacy of targeted therapies in cancer [36] and they can 
potentially address intra-patient heterogeneity. It wUl be necessary 
to pursue further studies not only of efficacy but also of toxicity 
before their use in patients can be considered. In this study, we 
have included human mononuclear cells and CD34+ HSCs as a 
part of toxicity studies (Table S7). In future studies it might be 
desirable to further optimize drug combinations against a range of 
control cells to ensure low toxicity of the drug combinations. Our 
approach also holds promise for the identification of personalized 
therapies, some of them composed of only two dugs, as shown in 
Figure 4 in the case of the PD325901 and imatinib combination in 
one if the BCR-ABL+ patients. 

Another challenge originates from the low predictive value of in 
vitro studies, even if performed with patient derived primary cells. 
[37] We have recently analyzed [38] the metabolomics of the 
human ALL bone marrow microenvironment. This might 
contribute, together with additional molecular studies addressing 
other components of this niche, for example cytokines, to the 
design of cell culture media that mimic more closely the patient 



conditions and accelerate the translation of our combinatorial 
therapy findings. 

In summary, using an experimental search algorithm we have 
identified drug combinations including imatinib. These drug 
combinations are effective in a cell line and in primary cells and 
may provide a first-step toward finding combination therapies for 
the treatment of leukemic patients, especially if BCR-ABL-H. 

Methods 

Ethical Statement 

All clinical investigations were conducted according to Decla- 
ration of Helsinki principles. All human studies were approved by 
the UCSD Human Research Protections Programs IRB. Written 
informed consent was received from participants prior to inclusion 
in the study. Written informed consent and written parental 
permission were obtained in accordance with Institutional Review 
Board guidelines. 

Cell lines and patient cells 

K562 and IMR-90 cells were maintained in RPMI (Invitrogen, 
USA) supplemented with 10% fetal bovine serum (FBS, Hyclone, 
USA), 100 U/ml penicillin, 100 |Xg/ml streptomycin (Hyclone, 
USA) at 37°C in 5% CO2. CD34-H hematopoietic stem cells were 
purchased from AllceUs and maintained in StemSpan SFEM 
(STEMCELL Technologies, BC, Canada) with 100 ng of SCR, 
TPO, G-CSF (Life Technologies, CA, USA), and 1 |xM of 
StemRegeninl (EMD millipore, CA, USA). 

Bone marrow specimens from ALL BCR-ABL- patients were 
obtained at the time of diagnosis from the Rady Children's 
hospital (San Diego, CA, USA). BCR-ABL-H ALL patient 
xenografts were obtained, as described in Bicocca et at., from the 
Muschen lab. [39] Diagnostic and cytogenetic information of the 
patients is included in Table S8. Blood from normal subjects was 
obtained from the Scripps Research Institute (San Diego, CA, 



Table 3. Results of statistical analysis using repeated measures one-way ANOVA with Dunnett's multiple comparisons between 
patient cells and normal mononuclear blood cells (Nl) on viability of the best drug combinations originating from the cell lines 
search and from pairs of single agents plus imatininb. 





Dunnett's multiple comparisons 


Mean Diff. (combinations) 


Adjusted P Value 
(combinations) 


Mean Diff. 
(pairs) 


Adjusted P 
Value (pairs) 


Nl vs. BOOl 


0.3468 


<0.0001 


-0.1379 


0.0427 


Nl vs. B002 


0.4064 


<0.0001 


-0.195 


0.077 


Nl vs. B003 


0.4570 


<0.0001 


0.6418 


0.9174 


Nl vs. B004 


0.3884 


<0.0001 


-0.8320 


<0.0001 


Nl vs. B031 


0.1914 


<0.0001 


-0.2496 


<0.0001 


Nl vs. B032 


0.314 


<0.0001 


-0.304 


<0.0001 


Nl vs. B033 


0.3054 


<0.0001 


-0.08711 


0.121 


Nl vs. B035 


0.3132 


<0.0001 


-0.03637 


0.951 


Nl vs. B036 


0.06369 


0.2464 


-0.1242 


0.0043 


Nl vs. B037 


0.258 


<0.0001 


-0.05527 


0.1199 


N1 vs. Burkitt 


0.1959 


<0.0001 


-0.03135 


0.7161 



Patients BOOl and B002 are BCR-ABL+ (ALL). Patients B031 to B037 are BCR-ABL- (ALL). 
doi:1 0.1 371/journal.pone.01 02221 .1003 
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USA). The study was performed in accordance with Institutional 
Review Board guidelines after obtaining informed consents from 
the patients. Briefly, mononuclear cells were isolated by FicoU- 
paque (GE Healthcare, CA, USA) density gradient centrifugation 
(1.077 g cm" ') at 400 g for 30 min without brake followed by 
three washes in PBS. Patient xenografts and ALL patient 
mononuclear cells were subsequendy cultured in RPMI with 
10% FBS. 

Kinase inhibitor library screen 

We identified candidate hits from the EMD kinase library 
(EMD MiUipore, USA) screening at 1 (xM concentration of kinase 
inhibitors with/without imatinib (0.125 \iM). The values of 
relative cell survival from four replicates of two independent 
screens were averaged. Selectivity (Vn/v„ where v,, and v^ indicate 
the viability of IMR-90 and K562) and imatinib combination 
index (v/vimatmib> whcrc v and Vimatimb indicate the viability of 
K562 without/with imatinib) were calculated and rank-ordered 
from the highest to the lowest. 

Combinatorial studies 

The combinatorial drug studies were done with the Echo 
automated liquid handler (Labcyte, CA, USA). Using our lab- 
developed high throughput screening manager software, the input 
files for combinatorial drug transfers were generated and drug 
droplets were transferred from the drug source plate to the 
designated wells of destination plates in 2.5 nl increments. 

Measurement of cell viability 

Cell survival was assessed by luciferase-based assay, ATPlite 
(PerkinElmer, CA, USA), which determines viable cell numbers by 
measuring the presence of ATP in all metabolicaUy active cells. 
For the measurement of cell viability, K562 and IMR-90 cells 
were plated at a density of 2000 and 1500 cells/ 384-well, 
respectively whereas primary patient cells were seeded at a density 
of 8000 cells/384-2311. Subsequently, the cells were treated with 
the drugs and 72 hours later, the ATPlite assay was performed 
according to the manufacturer's protocol, and luminescence was 
read with an Analyst HT instrument (Molecular Devices, CA, 
USA). IC50 values were calculated from nonlinear regression 
analysis using Prism 6 (Graphpad). 

Gene expression correlation analysis 

The Pearson correlation and p-values between mRNA expres- 
sion and drug IC50 values for each gene were computed over 789 
cell lines for 132 drugs [15]. This computation was repeated for all 
genes identified on the microarray. Next, the overlap of significant 
genes (p-value <0.001) between drugs and imatinib was computed 
and separated into four distinct comparisons. Given that four 
combinations of two drugs with either positive or negative 
correlations (namely negative-negative, negative-positive, posi- 
tive-negative, and positive-positive c;orrelation pairs) for a given 
gene are possible, we computed the number of shared genes for 
each possible drug pair and each of the four sign pairings. Each 
gene could receive only one correlation and corresponding sign 
(positive or negative) for a given drug. Since some drugs have a 
larger number of significant genes, we re-weighted the overlap 
parameter in order to get a normalized value. The formula 

Overlapcount / min(sig(drugl ) ,sig(drug2)) 

was used to re-weight the overlap between gene-to-drug combi- 
nations, where overlap is the number of genes meeting the 0.001 



p-value significance cutoff that appear with desirc'd sign pair (e.g. 
pos-pos) for the two drugs tested. The denominator is the 
minimum number of significant genes meeting the p-value cutoff 
for either of the two drugs tested individually. 

Initial conditions for the search algorithms 

The initial set of multidrug combinations tested in the first 
iteration of the algorithm was obtained using results of the single 
agent and pair screening. Imatinib was included in all combina- 
tions at a fixed dose (0.125 |iM) and we considered all the possible 
pairs of the other eleven drugs of Table 1 at the doses '0' to '4' 
corresponding to concentrations given in Table 2 (a). The 
selectivity of a multidrug combination can be described as the 
function 

where N= 11 and rf, is the dose of the drug i. We can then estimate 

the selectivity function as a geometric average of the measure- 
ments in the single and pair screening as 

s(dududu-,dM)- y n,Sf-2(4) 

were Sij(di,di) indicates the selectivity measured using only drugs i 
and J and S, (fi?, )is the selectivity measured using only drug i. Note 
that since d can take the value '0', these expressions include also 
the selectivity values measured with Imatinib alone and Imatinib 
with a second drug. 

The estimates for all possible combinations based on this 
geometric average were then ranked and we used the top 10 
combinations as initial position of the "particles" in the "particle 
swarm optimization" described below. 

Algorithmic method 

The algorithm implemented in this search integrated a particle 
swarm optimization (PSO) method and a genetic algorithm (GA) 
[18]. We can think about the PSO in terms of a set of particles that 
move in the space of drugs and doses. The algorithm depends on 
two parameters, a and (3 that provide the relative contribution of 
the local and the global maximum to the velocity of the particles. 
We used 10 particles in our implementation and we set the 
maximum number of new combinations to be measured at each 
iteration to Max_comb = 270. The parameters of the algorithms 
were optimized through testing on simulated selectivity sets 
generated combining data from many previous multi-drug 
screenings on cell lines. The procedures for the PSO and GA 
implemented in the search are schematically described in BOX 1 
and BOX 2. (Figure S2) 

Statistical analysis 

Statistical analyses were performed using Prism 6 (Graphpad). 
Matched patient and normal subject specimens were compared 
using Dunnett's multiple comparison tests after repeated-measures 
analysis of variance (ANOVA). 

Supporting Information 

Figure SI 

(PDF) 

Figure S2 

(TIF) 
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Table SI Screen of EMD kinase inhibitor library. Kinase 
inhibitors were rank-ordered according to the imatinib combina- 
tion index (viability of K562/viability of K562 with imatinib). Top 
three kinase inhibitors were selected accordingly. 

(XLS) 

Table S2 Correlation sign analysis between imatinib 
and drug responses. Shared genes from positive /positive 
Pearson correlation between mRNA expression, imatinib IC50 
values and mRNA expression, drug IC50 values. 

(XLSX) 

Table S3 IC50 values of the selected kinase inhibitors. 

IC50 values ((xM) were calculated from the nonlinear regression 

analysis. 

PCLS) 

Table S4 Rank ordered selectivity when algorithm- 
generated combinations were tested at each iteration. 

PCLS) 

Table S5 The list of top 49 combinations and corre- 
sponding single doses. 49 drug combinations were selected 
using two criteria: top combinations with a 70% viability of IMR- 
90 and low toxicity combinations with a 90% viability of IMR-90. 

(XLSX) 
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